Observation of the Y(l 3 Dj) bottomonium state through decays to 7r + 7r Y(1S) 
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Based on 122 x 10 6 Y(3S) events collected with the BABAR detector, we have observed the T(l 3 Dj) 
bottomonium state through the Y(3S) —¥ 77T(1 3 Dj) —t r y r yir + ir~Y(lS) decay chain. The signif- 
icance for the J — 2 member of the X(l 3 Dj) triplet is 5.8 standard deviations including sys- 
tematic uncertainties. The mass of the J — 2 state is determined to be 10164.5 ± 0.8 (stat.) ± 
0.5(syst.) MeV/c 2 . We use the tv + ty invariant mass distribution to confirm the consistency of the 
observed state with the orbital angular momentum assignment of the Y(l 3 Dj). 

PACS numbers: 13.25.Hw, 14.40. Nd, 14.65.Fy 



Heavy quark bound states below open flavor thresh- 
olds provide a key probe of the interactions between 
quarks. The mass spectrum and branching fractions of 
these states can be described by potential models and 
quantum chromodynamics S-wave and P-wave 

bottomonium (bb) states were first observed in the 1970s 
and 1980s. Only recently Q| has a D-wave bottomo- 
nium state, the triplet Y(l 3 Dj) been observed, where 
J = 1,2,3. The separation between the members of the 
triplet (intrinsic widths about 30 keV/c 2 ) is expected to 
be on the order of 10 MeV/c 2 |2j. A single state, inter- 
preted to be the J = 2 member of the T(l 3 Dj) triplet, 
was observed Q by the CLEO Collaboration in the radia- 
tive T(l i D2) — > jjT(1S) decay channel, but the quan- 
tum numbers L, J j|| and parity P were not verified. 

In this paper, we report the observation of the J = 2 
state of the T(l 3 Dj) in the hadronic w + ir~Y(lS) de- 
cay channel, with Y(1S) -)■ £+i~ (£ = e, //). This 
decay channel has been of interest for decades 0, @- 3 ■ 
Predictions for the branching fraction vary widely 6j- 
8]. It provides better mass resolution than the -fjY(lS) 
channel and allows L, J, and P, for which there is 
currently no experimental information, to be tested, 
through measurement of the angular distributions of the 
tt^~ and The only previous result for this chan- 

nel is the 90% confidence level (CL) branching fraction 
upper limit Bt^s^-iT^d ,) x B T ^3 Dj ^+„- r ^ S ) X 



B 



- < 6.6 x 10" 6 



r(is)-ye+t 

The analysis is based on a sample of (121.8 ± 1.2) x 10 6 
Y(3S) decays collected with the BABAR detector at the 
PEP-II asymmetric-energy e + e~ storage rings at the 
SLAC National Accelerator Laboratory, corresponding 
to an integrated luminosity of 28.6 fb . The BABAR 
detector is described elsewhere [9]. Monte Carlo (MC) 
event samples that include simulation of the detector re- 
sponse are used to determine the signal and background 
characteristics, optimize selection criteria, and evaluate 
efficiencies. Pure electric-dipole transitions [13] are as- 
sumed when generating radiative decays. 

The Y(l 3 Dj) in our study are produced through 



Y(3S) 7Xfcj'(2P) -> 77r(l 3 P>j) transitions, with 
J' = 0,1,2. To reconstruct the Y(3S) -> jjir+ir- £+£- 
final states, we require exactly four charged tracks in an 
event, two of which are identified as pions with oppo- 
site charge and the other two as either an e + e~ or /x + /x~ 
pair. Pion candidates must not be identified as electrons. 
To reject Bhabha events with bremsstrahlung followed 
by 7 conversions, we require the cosine of the polar an- 
gle of the electron with respect to the e~ beam direc- 
tion to satisfy cos0 e - < 0.8 in the laboratory frame. To 
improve the energy measurements, up to three pho- 
tons are combined with e candidates to partially recover 
bremsstrahlung [ll|. The Y(1S) candidate is selected 
by requiring —0.35 < m e + e - — ?7ir(is) < 0.2 GeV/c 2 or 
\ m ti+^- ~ m r(is)\ < 0-2 GeV/c 2 , where the invariant 
mass of the lepton pair rrip+ p- i s then constrained to the 
nominal Y(1S) mass value [1 21 ] . The pion pair is com- 
bined with the Y(IS) candidate to form a Y(l 3 Dj) can- 
didate (mass resolution 3 MeV/c 2 ). To eliminate back- 
ground from 7 — > e + e~ conversions in which both the e + 
and e~ are misidentified as pions, we reject events with a 
cosine for the laboratory 7r + 7r _ opening angle cos#„.+ i7r - 
greater than 0.95 if the converted e + e~ mass is less than 
50 MeV/c 2 , and events with a laboratory angle between 
the 7r + 7r~ pair and £ ± that satisfies cos 6 ■ K + 7T - g± > 0.98. 

Photons from Y(3S) -> 7 XfcJ'(2P) (xw(2P) -> 

7T(l 3 Dj)) decays have energies between 86 and 

122 MeV (80 and 117 MeV @|) in the Y(3S) center- 
of-mass (CM) frame. Our resolution for 80 MeV pho- 
tons is about 6.6 MeV. We require at least two pho- 
tons in an event: one (the other) with CM energy larger 
than 70 MeV (60 MeV). Photons from final-state ra- 
diation (FSR) are rejected by requiring the cosines of 
the laboratory angles between the cascade photons and 
leptons to satisfy cos#£ i7 < 0.98. In case of multi- 
ple photon combinations, we choose the one that min- 



imizes x 2 — ^2 {e. 

=1,2 



El 



'/o-| f , where E^ p 



are the 

7X&j'(2P)) or expected 3 



7 "oxp; / " /<; 
nominal [ij (for Y(3S) -> 
(for x&J'(2P) 7^(1 3 -D./)) photon energies that corre- 
spond to one of the six possible Y(3S) — » 7XbJ' (2P) — I 



j-fT(l 3 Dj) transition paths allowed by angular momen- 
tum conservation, with [a E i ) the measured energies 
(resolutions). We verified that the x 2 procedure does not 
bias our results, using simulated data samples in which 
the assumed Y(l 3 Dj) mass values are varied. 

The T(l 3 Dj) candidate is combined with the two pho- 
tons to form a T(3S) candidate, whose CM momentum 
must be less than 0.3 GeV/c. The T(3S) mass is then 
constrained to its nominal value [ljj]. The Y(3S) labora- 
tory energy (resolution 25 MeV) is required to equal the 
summed e + and e~ beam energies to within 0.1 GeV. 

We identify four background categories within our fit 
interval 10.11 < m^- l+l - < 10.28 GeV/c 2 : Y(3S) 
decays to (I) 7X6./' (2P) with XbJ'(^P) ->• uiT(1S) and 
uj -> 7r+7r-(7r°), (II) ir+ir-T(lS) with FSR, (III) r)Y(lS) 
with n -> 7r+7r-7r°(7), and (IV) ^T{2S) or 7r°7r°r(2S) 
with T(25) — > ir + n~T(lS). Categories I and II are the 
main backgrounds. 

An extended unbinned maximum likelihood fit is ap- 
plied to the sample of 263 events in the fit interval. The 
fit has a component for each of the three T(l 3 Dj) signal 
states and four background categories. The likelihood 
function is C — exp (— X) 3 - rij ) Ilili V , " ,Pj ■in,) , 
with N the number of events, rij the yield of compo- 
nent j, Vj the probability density function (PDF) for 
component j, and m the n + Tr~ £ + £~ invariant mass. 

The PDFs are derived from MC simulations. Each 
signal PDF is parameterized by the sum of two Gaus- 
sians and a Crystal Ball (CB) function For back- 

ground category I, we use the sum of a CB function, 
which describes the lu — > 7r + 7r~7r° events, and two Gaus- 
sians, which model the two peaks from Xbi,2(2P) decays 
to uT(lS) with uj — > tt + tt~. A bifurcated Gaussian, 
a high statistics histogram, and a Gaussian, model the 
PDFs for backgrounds II, III, and IV, respectively. 

A large data control sample of T(3S) — > jXbJ'i^P) ~> 
JjT(2S) events with T{2S) -> 7r+7r-T(lS') and T(IS) -> 
£ + £~ is used to validate the signal PDFs and mass recon- 
struction. The control sample is selected using similar 
criteria to those used to select the T(l 3 Dj). The back- 
ground contamination is about 2%. Only a small dif- 
ference is observed between the shapes of the T(2S) —> 
7r + 7r~ £ + l~ invariant mass distributions in data and simu- 
lation. The signal PDF is adjusted to account for this dif- 
ference. The reconstructed T(2S) mass is shifted down- 
wards by 0.70±0.15 (stat.) MeV/c 2 compared to its nom- 
inal value [12J . We apply this shift as a correction to the 
T(l 3 Dj) mass results presented below. 
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FIG. 1: The ir + n~£ + £~ mass spectrum and fit results. The 
two peaks near 10.25 GeV/c 2 arise from Xbj>(2P) — > uY(lS) 
background events with cj — > ir + n~ . 



Eleven parameters are determined in the fit: the three 
signal yields and three masses, the yields of background 
categories I and II, and - within background category I - 
the Xbi{2P) mass and the relative yields of the Xfci(2-P) 
and Xb2(2P) peaks from u — > it + tt~ decays. The mass 
difference between the Xbi (2P) and Xb2{2P) peaks is 
fixed to its measured value [12j . The yields of background 
categories III and IV are fixed to their expected values 

13 IS 1 . 



based on the measured branching fractions 

Figure Q] shows the ir + ir~£ + £~ mass distribution and 
fit results. The results for the separated T(1S) — > e + e~ 
and T(1S) — > n + fJ.~ channels are shown in Fig. [2] The 
e + e~ channel has a smaller efficiency than the 
channel in part because of the criteria to reject Bhabha 
events. The differences in efficiency between the e + e~ 
and channels, including those for the XbJ'(2P) —> 

ujT(1S) background events, are consistent with the ex- 
pectations from the simulation within the uncertainties. 
We find 10.6±|;| T^Lq), 33.91^ T(l 3 £> 2 ), and 9.4^; 2 
T(1 3 D3) events. The positions of the three signal peaks 
in Fig. [T] are stable with respect to different initial as- 
sumptions about their masses within the fit interval. The 
fluctuations at around 10.13 and 10.18 GeV/c 2 are dis- 
cussed below. The fitted background category I and II 
yields of 50 ± 9 and 94 ± 13 events agree with the MC 
expectations of 51 and 94 events, respectively. The fitted 
Xbi(2P) mass value of 10255. 7±0.7 (stat.) MeV/c 2 (after 
applying the shift of +0.7 MeV/c 2 from the T(2S) mass 
calibration) is in good agreement with the nominal value 
10255.5 ± 0.5 MeV/c 2 12], validating the calibration. 

Fit biases are evaluated by applying the fit to an en- 
semble of 2000 simulated experiments constructed by 
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10.15 



10.2 



10.25 



jt+ji"e + e" mass (GeV/c ) 




10.15 



10.2 



10.25 



Jt+Jiu+u: mass (GeV/c ) 



FIG. 2: The -k + tt £ + £ mass spectra for the separated (a) 
Y(1S) e + e~ and (b) Y(1S) -> channels. The results 

of the fit are shown. The legend is given in Fig. [1] 



randomly extracting events from MC samples. The num- 
bers of signal and background events and the T(l 3 Dj) 
masses correspond to those of the fit. The biases are 
1.6±0.1, -1.8±0.2, and 1.0±0.1 events for the T(l 3 L>i), 
Y(1 3 D 2 ), and Y(1 3 D 3 ), respectively. We subtract these 
biases from the signal yields. The biases on the masses 
are negligible. 

Multiplicative systematic uncertainties arise from the 
uncertainty in the number Nr(3S) of T(3S) events in the 
initial sample (1.0%) and in the reconstruction efficien- 
cies for tracks (1.4%), photons (3.0%), and particle iden- 
tification (2.0%). Additive systematic uncertainties orig- 
inate from signal and background PDFs, evaluated by 
varying the PDF parameters within their uncertainties, 
background yields, evaluated by varying the background 
category IV (III) yield by its uncertainties (by ±100%), 
the fit bias, and the T(2S) mass calibration. The fit bias 
uncertainties are defined as the quadratic sum of half the 
biases and their statistical uncertainties. The mass cal- 
ibration uncertainty is taken to be half the Y(2S) mass 
shift added in quadrature with the Y(2S) mass uncer- 
tainty 12]. The overall additive uncertainties for the sig- 



nal yields (masses) are 1.5 — 2.0 events (0.48 MeV/c 2 ) and 
are dominated by the contribution from the background 
yields (Y(2S) mass calibration). 

As a check, we repeat the fit with an additional back- 
ground term, given by a second order polynomial. The 
purpose of this check is to test for the effect of potential 
unmodeled background. The parameters of the polyno- 
mial are left free in the fit (thus there are 14 free pa- 
rameters). The fitted Y(l 3 Dj) yields are affected by less 
than 0.5 events compared to our standard fit, for all J 
values. The shifts in the fitted mass values are less than 
0.05 MeV. Since this polynomial background term is not 
motivated by any known source and since the description 
of the background without the additional term is good, 
we do not use this alternate background model to define 
a systematic uncertainty. 

We define the statistical significance of each Y{l 3 Dj) 
state by the square root of the difference between the 
value of — 21n£ for zero signal events and assuming the 
bias-corrected signal yield, with the masses and yields of 
the other two states held at their fitted values. These 
results are validated with frequentist techniques. Sys- 
tematics are included by convoluting C with a Gaussian 
whose standard deviation (cr) equals the total systematic 
uncertainty. The significances of the T(l 3 Di), X(1 3 Z?2), 
and Y(1 3 D 3 ) observations are 2.0 (1.8), 6.5 (5.8), and 
1.7 (1.6) a without (with) systcmatics included, respec- 
tively. If we use the raw signal yields, rather than the 
bias-corrected yields, the statistical significances of the 
J = 1, 2 and 3 states are 2.4, 6.2, and 2.0 cr, respectively. 

From Fig. [T] it is seen that the data exhibit up- 
ward fluctuations at ir + n~£ + £~ masses around 10.13 and 
10.18 GeV/c 2 . To investigate the significance of these 
fluctuations, we re-perform the fit with the J = 1 mass 
constrained to 10.13 GeV/c 2 rather than leaving it as a 
free parameter. An analogous fit is made with the J = 3 
mass constrained to 10.18 GeV/c 2 . The statistical signif- 
icance for this alternate J = 1 (J = 3) peak, evaluated 
using the raw signal yield, is 2.0 cr (1.3 cr), compared to 
2.4 a (2.0 cr) for our standard fit. The J — 2 signal yield 
and mass shift by less than 1 event and 0.04 MeV/c 2 , 
respectively, in these alternate fits. 

We determine branching fractions by dividing the bias- 
corrected signal yields by the selection efficiencies and 
Ntqs). The significances of the T(l 3 Di) and Y(1 3 D 3 ) 
peaks are low and we do not have clear evidence for 
them. For the J = 1 and 3 states, we also present up- 
per limits on the branching fractions assuming the fitted 
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FIG. 3: (a) The tt + tv~ mass spectrum in the T(1 3 D2) sig- 
nal region. The area under each curve equals the number of 
events. (b,c) Distributions in the T(1 3 Z>2) signal region of (b) 
the angle x between the tv + tt~ and £ + £~ planes, and (c) the 
7r + helicity angle. The uncertainties include both statistical 
and systematic terms. 



masses. The efficiencies for the six allowed T(3S) — > 
7Xbj'(2P) ->■ 77T(1 3 A 7 ) paths differ by up to 7.5% and 
therefore do not factorize, leaving six unknown branching 
fractions but only three measured signal yields. How- 
ever, 91.4% of the T{3S) -)• 77 :T(l 3 L>i) and 88.7% of 
the T(3S) — > 77T(1 3 Z?2) transitions are predicted [1] 
to proceed through the Xbi(2P) state, while T(3S) — > 
77T(1 3 Z)3) transitions can only proceed through the 
Xf,2(2P). Therefore, we evaluate the branching fractions 
for the dominant modes only, using the predicted ratios of 
the branching fractions to account for the non-dominant 
transitions. The efficiencies of the dominant modes, av- 
eraged over the T(1S) — > e + e~ and final states, 
are 26.7 ± 0.1%, 26.7 ± 0.1%, and 25.7 ± 0.2% for the 
T(l 3 L>i), T(1 3 L> 2 ), and T{1 3 D 3 ), respectively. 

The branching fraction products for the dominant 
modes Bjij = Br(3S)^x bJ '(2P) x B Xbj ^ 2 p)^r(i 3 Dj) X 



B 



x B 



(or the upper limits at 



90% CL with systematics included) are, in units of 10 7 , 
B n = 1.27±g;|| ± 0.28 (< 2.50), B 12 = 4.91^ ± 0.3, 
and B 23 = 1.34±g |l ± 0.24 (< 2.80). We determine the 
T(1 3 D 2 ) mass to be 10164.5 ± 0.8 ± 0.5 MeV/c 2 , which 
is consistent with, and more precise than, the result 
10 161.1 ± 0.6 (stat.) ± 1.6 (syst.) MeV/c 2 from CLEO 

From the T(3S) — > 'JXbJ' (2-P) branching fractions and 
uncertainties [12j and XM'i^P) ~~ lY(l 3 Dj) branch- 
ing fraction predictions [I] we determine B[T{l 3 Dj) -» 
n + Tr~Y(lS)] (or 90% CL upper limits including system- 
atics) to be 0.42 ±£^±0.10% (< 0.82%) for the T^Di), 
0.66 ±£^±0.06% for the T{1 3 D 2 ), and 0.29 ±£; 22 ±0.06% 
(< 0.62%) for the T(1 3 D 3 ), which lie between the predic- 
tions of about 0.2% from Ref. Q and 2% from Ref. Q. 

Figure [3Ja) shows the 7r + 7r~ mass distribution 



for events in the T(l 3 £>2) signal region 10.155 < 
rrin+Tr-i+e- < 10.168 GeV/c 2 after subtraction of the 
backgrounds using the estimates from the fit. The data 
are corrected for mass-dependent efficiency variations. 
Shown in comparison are the expectations for the de- 
cay of a D 15j, 5* 15|, or 1 P\ 16[ bottomonium state to 
tt+tt-T(1S). The resulting X 2 probabilities of 81%, 11%, 
and 10%, respectively, strongly favor the D state. 

The distribution of the angle \ between the £ + £~ and 
7r + 7r~ planes in the T(l 3 Dj) rest frame, for events in the 
T(1 3 D2) signal region, is shown in Fig. EJb). The data 
are corrected for background and efficiency. The x dis- 
tribution is expected to have the form 1 + j3 cos 2% with 
sign(f3) = (-1) J P [ill], where P is the parity. A fit to the 
data yields (i = -0.41 ± 0.29 (stat.) ± 0.10 (syst.), consis- 
tent with the expected assignments J — 2 and P = — 1. 

The background-subtracted, efficiency-corrected dis- 
tribution of the helicity angle 9^, for events in the 
T(1 3 D2) signal region, is shown in Fig. [3jc), where 0^ 
is the angle of the 7r + in the tt + tt~ rest frame with re- 
spect to the boost from the T(1 3 D2) frame. For D-state 
decays to ir + ir~ -T(1S), 6 n follows a 1 + |(3 cos 2 n - 1) 
distribution, where £ is a dynamical parameter to be 
determined experimentally. For S'-state decays, the 9^ 
distribution is flat (£ = 0). A fit to data yields £ = 
— 1.0 ± 0.4 (stat.) ± 0.1 (syst.), disfavoring the S state. 

In summary, we have observed the T(1 3 I?2) bottomo- 
nium state through decays to 7r + 7r~X(lS'). The signif- 
icance is 5.8cr including systematic uncertainties. We 
improve the measurement of the T(l 3 £)2) mass and de- 



termine the T(l 3 Dj) 



•7T T(1S) branching fractions 



or set upper limits. We use the 7r + 7r~ invariant mass, 
the angle between the 7r + 7r~ and £ + £~ planes, and the 
7r + helicity angle, to test the consistency of the observed 
state with the expected quantum numbers L = 2. J = 2 
and parity P = — 1 for the dominant member of the 
triplet T(1 3 D 2 )- 
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